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ABSTRACT 

Context. The basic properties of galaxies can be affected by both nature (internal processes) or nurture (interactions and effects of 
environment). Deconvolving the two effects is an important current effort in astrophysics. Observed properties of a sample of isolated 
galaxies should be largely the result of internal (natural) evolution. It follows that nurture-induced galaxy evolution can only be 
understood through comparitive study of galaxies in different environments. 

Aims. We take a first look at SDSS (g-r) colors of galaxies in the AMIGA sample involving many of the most isolated galaxies in the 
local Universe. This leads us to simultaneously consider the pitfalls of using automated SDSS colors. 

Methods. We focus on median values for the principal morphological subtypes found in the AMIGA sample (E/S0 and Sb-Sc) and 
compare them with equivalent measures obtained for galaxies in denser environments. 

Results. We find a weak tendency for AMIGA spiral galaxies to be redder than objects in close pairs. We find no clear difference 
when we compare with galaxies in other (e.g. group) environments. However, the (g-r) color of isolated galaxies shows a Gaussian 
distribution as might be expected assuming nurture-free evolution. We find a smaller median absolute deviation in colors for isolated 
galaxies compared to both wide and close pairs. The majority of the deviation on median colors for spiral subtypes is caused by a 
color-luminosity correlation. Surprisingly isolated and non-isolated early-type galaxies show similar (g-r). We see little evidence 
for a green valley in our sample with most spirals redder than (g-r)=0.7 having spurious colors. 

Conclusions. The redder colors of AMIGA spirals and lower color dispersions for AMIGA subtypes -compared with close pairs- is 
likely due to a more passive star formation in very isolated galaxies. 

Key words, galaxies: general — galaxies: fundamental parameters — galaxies: interactions — galaxies: evolution 



1. Introduction 

In some ways our study of galaxy properties is still in its infancy. 
We have a plethora of theoretical ideas but we are still trying to 
understand what basic measures e.g. size, luminosity and color 
are telling us. This is, in part, because galaxies are composite 
structures and in another part because we cannot easily separate 
effects of nature (internal processes) from nurture (interactions 
and effects of environment). 

The optical colors of galaxies reflect their stellar populations 
and these colors correlate with morphology and environment. 
The distribution of g alaxy colors in the (g-r) vs. (u-g) plane 
dStrateva et al.ll200"ll) shows a strong bimodality with clear sep- 
aration into red and blu e sequences. Study of morphology and 
spectral classification ( Strateva et al. 2001) for subsamples of 
287 red and 500 blue galaxies show that the two color peaks cor- 
respond roughly to early- (E, SO, and Sa) and late-type (Sb, Sc, 
and Irr) galaxies as expected from the respective dominance of 
old and young stellar populations. Colors of galaxies also depend 
strongly on luminosity in the sense that mor e luminous galaxie s 
of the same morphological type are redder (Bal drv et al. 12004). 
The color-magnitude relation is most obvious in the rest-frame 
(g-r) (corrected to z=0.1), wher e the separation betw een the red 
and blue populations is evident (Blanton et al. 20031). 

Environment is also thought to play a role in the mix of 
morphological types for a sample of galaxies which is re- 
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fleeted by the morpho logy-density relation dDressleti [l980; 
Ivan der Wei et alJl2010l and references therein). In dense envi- 
ronments luminous red early-type galaxies predominate while in 
lowest density e nvironment blue late-type spirals a re the defin- 
ing population (lDressletHl980l ICapak et alll2007l) . While it is 
easy to recognize a rich cluster, definitions of low density envi- 
ronments can be confusing. In recent years there has been an in- 
creased emphasis on identifying low density or isolated galaxy 
populations. One of the most useful samples remains the visu- 
ally selected Catalo g of Isolated Galaxies (CIG) compiled by 
Karache ntseyal (1 1 9731) more recently vetted as the AMIGA sam- 
ple dSulentid l2010, and references therein). 

Star formation is strongly dependent on environ ment with in- 
creased activity towards low density environmen ts dLewis et al.l 
| 2002tlHogget al]l20 03. 2004: lBaldrv et al.ll2006l) . In this sense, 
Baloghetajj (1998) found that the mean star formation rate in 
galaxies with similar bulge-to-total (B/T) luminosit y ratios is 
alway s lower in clusters than in the field. However Patton et al. 
(1201 ll) found that the opposite is happening in galaxy pairs with 
clear signs of star formation induced by interaction within the 
blue galaxies. This result a grees with the bluer (U- B and B-V) 
colors found previously bv lLarson & Tinslevl d 19781) for peculiar 
galaxies. 

Approximately 10% of the low density universe is populated 
by pairs (also by triplets and dense groups) where environmen- 
tal effects can reproduce t he processes that happen in starburs t 
galaxies and clusters (e.g. IXu & Sulentid [l99lt Ilovinoll2002l) . 
What we are trying to say is that there are two sources of en- 
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vironmental effects: 1) the morphology-density relation and 2) 
one-on-one galaxy interactions which can be found in a wide 
range of density environments. In this sense, the AMIGA vet- 
ting of the CIG (see below) is intended to produce an isolated 
sample that minimizes both effects. 

The AMIGA project (Analysis of the interstellar Medium 
of Isolated GAlaxies) is producing and analyzing a mul- 
tiwavelength database f or a refinement of the Catalog of 
Isolated Galaxies (CIG, iKarachentseval IT973L n=1050 galax- 
ies). We have evaluated and improved the sample in dif- 
ferent ways: 1) Revision of all CIG positions in the sky 
dLeon & Verdes-Monten egro 2003), 2) optical characterization 
of the sample inc luding completeness, luminositi es, and helio- 
centric velocities ( Verdes-Mo ntenegro et al.ll2"005l and section 2 
of this paper), 3) Morpholo gical refinement and identification of 
galaxies with asymmetries dSulentic et alj|200q, and section 2.2 
of this paper), and 4) quantification of the degree of isolation 
including the number density to the 5th neighbour ( 77^) and the 
estimation of the tidal force (Q) dVerlev et alj|2007alfbb . 

The multiwavelength analysis of our AMIGA sample has 
shown that isolated galaxies have different properties than 
galaxies in higher density environments (even field samples). 
Variables expected to be enhanced by interactions are lower in 
AMIGA than any other sample, such as lower in frared luminos- 
ity (L fir< 10.5 Lo) and colder dust temperature (Lisenfel d et alJ 
2007), a low level of radio continuum emiss ion dominated 
by mild disk star formation dLeon et alJ 12008). no radio ac- 
tive galactic nuclei (AGN) selec ted using the radio-FIR cor- 
relation (0%; Sabater et al. 2008) and a small number of opti- 
cal AGN (22%: ISabater et alJl20l1 su bmitted), a smaller frac- 
tion of asymmetric HI profiles (< 20%. | Espada et al.ll201lh . and 
less molecular gas dLisenfeld et al.1 1201 lb . In addition, early- 
type galaxies in AMIGA are fainter than late types, and more 
AMI GA spirals hos t pseudo-bulges rather than classical bulges 
dDurbala et al.l 2008). The data are being released and periodi- 
cally updated at http://amiga.iaa.es, where a Virtual Observatory 
compliant web interface with different query modes has been 
implemented. 

This paper presents a first look at colors for the AMIGA sam- 
ple. It now becomes possible to analyze more than half of the 
AMIGA sample using uniform digital images, magnitude s an d 
colors of the sloan digital sky survey (SDSS. lYork et alj|2000h . 
In Sect. 2, the data revision of the AMIGA sample is provided. 
The sample selection is presented in Sect. 3, and the determi- 
nation of the absolute magnitudes is described in Sect. 4. The 
dependence of the rest-frame color as function of the morpho- 
logical type and environment is analyzed in Sect. 5 and 6, re- 
spectively. Finally, the conclusions are presented in Sect. 7. 

2. The AMIGA data revision 

Given that, since the project started, a significant number of new 
data are available for the sample, we have revised the appar- 
ent magnitudes, morphological types, distances and optical lu- 
minosities for the CIG sample w ith r espect to the ones us ed in 
IVerdes-Montenegro et ail d2005l) and lLisenfeld etalJ d2007l) . In 
the frame of the Wf4Ever project 1 , the revision of these proper- 
ties has been partially automated by the implementation of sci- 
entific workflows which ask for and gather values, of some of 
the properties for all the AMIGA galaxies, from the HyperLeda 

1 http://www.wf4ever-project.org 

2 http://www.taverna.org.uk/ 

3 http://www.myexperiment.org/packs/23 1 .html 



catalog (see details below). These workflows allow also com- 
parison between values stored in the AMIGA database coming 
from different releases, registering the new values while keeping 
a version of the old ones. Therefore the methodology to calculate 
these values is stored in these workflows enabling reproducibil- 
ity and re-usability, even the re-purposability for similar cross- 
boundary use cases. Currently releases of data in HyperLeda 
are not registered, hence while the methodology, as explained 
above, can be preserved, traceability of the HyperLeda data is 
not possible yet. These workflows have been built in the Taverna 
Workflow Mamagement System 2 and can be accessed in the 
My Experiment portal 3 . The data for the full AMIGA sample are 
listed in Table Q] and are also available in the AMIGA VO com- 
pliant interface. We describe the details of the revision below. 
For completeness, we explain here all the revised parameters, 
although optical luminosity has not been used in this paper. 

2.1. Velocities and Distances 

We have searched for new velocities for those CIG with helio- 
centric velocities (Vhei) larger than 1000 km/s, since for lower 
values redshift-independ ent distance estimates have been pre- 
ferred, and the same as in lVer des-Monteneg ro et al] (H005) were 
kept. The value of for a galaxy with available redshift data 
has been updated when the errors quoted in the recent bibliogra- 
phy were smaller than in our previous compilation. For a subset 
of galaxies we foun d that the source of t he data provided in NED 
were obtained from Stocke et al. I d2004h but, while the paper was 
providing galactocentric velocities, NED was quoting them as 
heliocentric. This error was reported to NED team and both cor- 
rected in NED and in our database. Only for 40 CIG galaxies 
velocity data were not found in the bibliography, hence 1010 out 
of 1050 CIG galaxies currently have redshift measures. 

Distances were obtained from virgo-centric velocities (V vl >) 
as D = Vvi r /Ho (using H ( > = 75 km/s Mpc~'), where V„> is cal- 
culated from the heliocentric velocity V/ K ; and galactic coordi- 
nates following HyperLeda convention. Virgo-centric velocities 
were calculated as Vhei = Vig + 208 cos(6*), where V^ is the 
radial velocity with respect to the Local Group and 9 is the an- 
gular distance between the galaxy and Virgo center. In order to 
transform the Vhei to V V i r we transf ormed a and 6 coordinates in 
iLeon & Verdes-Montenegrol d2003l) to I and b values. 

2.2. Morphologies 

ISulentic et al] d2006l) performed a careful revision of the mor- 
phologies for the whole CIG sam ple based on POS S II (Second 
Palomar Observatory Sky Survey; Rei d et alJI 19911) images. We 
have recently revised those morphologies for all CIG galaxies 
with Vhei > 1000 km/s and CCD images available either from 
SDSS or our own data (N = 843). For CIG galaxies with V he i < 
1000 km/s (N = 57), morphological types were compiled from 
the bibliography and the mean value used. For N = 134 galax- 
ies with only POSS II data available, a second revision has been 
performed, with mi nor modifications to the values assigned in 
ISulentic et alJ d2006l) . For the remaining galaxies for which we 
were unable to perform any classification, we used data found in 
NED and HyperLeda database. 

Globally there is a shift in the morphologies with respect to 
ISulentic et al.l (2006) towards later types by AT = 0.2, that we in- 
terpret as due to the higher resolution provided by CCD images 
revealing smaller bulges. Types (T) have b een coded fol l owing 
RC3 (see morphological codes in Tab. 0. S ulentic et al.l {2006) 
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Table 1. Main properties of the CIG sample 1 . 
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1 The full table is available in electronic form at http://amiga.iaa.es The columns correspond to (1) Galaxy identification according to CIG 
catalog; (2) Heliocentric velocity; (3) Distance; (4) Morphological type; (5): Degree of optical asymmetry, (6) B-band magnitude from 
HyperLeda, (7): Galactic extinction, (8): Internal extinction; (9) K correction in the B-band; (10): Corrected magnitude in the B-band; (11) 
logarithm of the optical luminosity in the B-band; (12): Major axis; (13) Minor axis; (14) Galaxy inclination, as described in Sect. 2. In 
column (15) we flag galaxies wich code 1 when the galaxy is part of the studied in this paper, and otherwise, following Sect. 3. 



also flagged galaxies suspected to be interacting. We have now 
replaced this code by a descriptive one, based also on visual in- 
spection of the optical images, as described next. Code IA=0 
is assigned if no relevant signs of distortions are visible. Code 
IA=1 corresponds to galaxies seen in the images as asymmet- 
ric, lopsided, warped or distorted, with an integral sign shape 
or tidal feature (tail, bridge, shell). Those galaxies are always 
treated separately in our analysis, in order to check whether they 
present a different behaviour than the ones classified with code 
0. Code IA=2 is ass igned if the galaxy was identified as pairs in 
IVerlev et al.l (l2007al) and/or in NED, or if the galaxy looks like a 
merger or superposition of two galaxies. These galaxies are not 
considered part of AMIGA sample and are marked differently in 
the plots. The new data are presented in Tab.Q]and now all CIG 
galaxies have a morphological code. 

2.2.1. Apparent magnitudes 

Blue optical m agnitudes B7- have been obtained from HyperLeda 
(Patu rel et al.l l2003). instead of using the Catalog of G alaxies 
and Clusters of Galaxies (|Zwickv et alJ I1961-1968I) as in 



IVerdes-Montenegro et al.l (12005 ). and the corrected was cal 
culated as follows: 



B^ — Bt 



A, - A* 



(1) 



where A g is the galactic dust extinction, A,- the internal ex- 
tinction correction and Ag the K correction. A g has been ob- 
tained from HyperLeda (ISchlegel et al.lll998l) . except for CIG 
447 for which only NED data were avail able. A; was calculated 
as in HyperLeda dBottinelli et al.l 0*9951) using their values for 
log(R25) (axis ratio of the isophote 25 mag/arcsec2 in the B- 
band for galaxies) but our own morphological types. Finally Ak 
was obtained as defined in HyperLeda 2 (A^=ak(T)V/ 1( ,;/10000; 
Ide Vaucouleurs eta i1 ll976l) . using our own morphologies and 
distances. Both uncorrected and corrected magnitudes are listed 
in Tab. □ 



2.3. Optical luminosity 

The optical luminosities have been derived as 

log(L B /L ) = 11.95 + 21og[D(Mpc)] - 0.4B^ 



(2) 



where the solar lumino sity is given in units of the solar bolo- 
metric luminosity as in LisenfeWetaD (120071) . The mean value 
of the difference between our n ew values for optical lu minosities 
with respect to the ones used in lLisenfeld et al.l ( 12007b is as small 
as 0.02 + 0.18 mag but has the advantage of being derived from 
Br magnitudes compiled from HyperLeda database, so reducing 
the errors. 



2.4. Optical diameter and axis ratio (D25 and R25) 

We have made use of the HyperLeda database in the compilation 
of the major axis D25 (isophotal level at 25 mag/arcsec 2 in the B- 
band) and axis ratios R25 {R25 = D25/1J25, where 1J25 is the minor 
axis, at an isophotal level of 25 mag/arcsec 2 ), as well as their 
errors (minor and major axis are given in Tab.[TJ. 



2.5. Inclinations 

The inclination (;') was determined from the value of R25 in 
HyperLeda database an d our revised morpho logical type with 
the following equation (iHeidmann et al.lll972b : 



skr(i) 



I _ jQ-21og(R 25 ) 

1 _ io- 21 °s(Ro) 



(3) 



http://leda.univ-lyon 1 .fr/leda/param/btc.html 



where log(fl ) = 0.43 + 0.053 x T for T< 7 and log(/? ) = 0.38 
for T > 7. 



3. Sample selection 

As starting sample we u s e the ca talog of 791 AMIGA galaxies 
selected by IVerlev et al.l d2007bl) . where the galaxies with iso- 
lation parameters Q>-2 and rj^lA (the tidal strength created 
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Fig. 1. Distribution of the heliocentric velocities of the 466 CIG 
galaxies used in this work. 



h of the internal bind- 
rjk, this translates to a 



by all the neighbours, Q, is more than 1 
ing forces; f or the local number density, 

value of 2.4; lAthanassoulall 1 9841) . and with recession velocities 
V r <1500 km/s, were rejected. These conditions imply that all 
galaxies in our starting sample h ave their evolution d ominated 
by their intrinsic properties (see Verl ev etal] |2007b. for more 
details). This sample is complete up to B^=15.3, where B C T is 
the magnitude in the B-band from Hy perLeda after corrections. 
The co mpleteness limit was decided in lVerdes-Montenegro et alJ 
(2005) (B^=15 mag) after applying the <V/Vm> test (iSchmidt 
1968), and has been recently updated since B-band magnitudes 
were recalculated following HyperLeda expressions in 2010 (see 
section 2). There are 657 objects in the complete AMIGA sam- 
ple that also fulfill the above isolation criteria. Photometric data 
used in this work come from the SDSS-III (Data Release 8, 
DR8 lAiharaet alJ|201l[) . The SDSS project used a 2.5 m tele- 
scope dGunn et al.l2 006) to obtain photometric information from 
CCD images in u, g, r, i, and z passbands. A new approach for 
background subtraction was applied in DR8 that first models the 
brightest galaxie s in each field in orde r not to affect the estimated 
sky background (Blanto n et alj|20l"il) . 

From the 657 galaxies that are part of the complete AMIGA 
sample, we found 496 objects in the SDSS database. We found 
and removed 14 galaxies from the catalog involving imprecise 
identifications of galaxy nuclei in the SDSS database or galaxies 
with a nearby bright star adversely affecting galaxy photome- 
try. We also removed 16 galaxies with unknown redshifts since 
those will be needed for the following analysis. The final sample 
consists of 466 isolated galaxies with radial velocities between 
1500 and 24000 km s -1 . Hereafter, we refer to it as the AMIGA- 
SDSS sample. We have represented its redshift distribution in 
Fig.Q] This sample includes 70 galaxies with asymmetries that 
are suspected of being involved in interactions or having nearby 
companions. These properties are represented in our sample by 
the degree of asymmetry IA=1 (63) or IA=2 (7) (see Section 
2.2), where IA=2 represents the most asymmetric galaxies. We 
do not reject these objects since we are trying to see what fea- 
tures can a galaxy develop in isolation, but, throughout the paper, 
we check the effect of these galaxies in the median colors. We 
focus our attention on the two most statistically significant sub- 
samples involving Sb-Sc spirals and E/S0 early-types. These in- 
volve approximately 68% and 15% of the AMIGA-SDSS sam- 
ple used in this work respectively. The former is interpreted as a 
kind of "parent population" of isolated galaxies because it rep- 
resents 2/3 of the total sample. The latter was somewhat of a 



surprise under the assumption that isolated galaxies are simply 
the outliers of loose grou ps where the early-type fraction is ~ 0.4 
( Ivan der Wei et alj|20 07). We find ~0. 15 in our sample, which is 
still large enough to be interesting. Other types are represented 
in such small numbers that we assume them to be simply ran- 
dom galaxies that happen to be unusually isolated at this time. 
This assumption is motivated by the standard idea that all disky 
galaxies with a large bulge are products of nurture. The late- 
type spiral populatio n involves galaxies with small bulges (or 
even pseudo-bulges, Durbalaet all l2008h which, if interpreted 
as a secular evolution diagnostic, cannot have spent much time 
in richer environments. 

4. Deriving rest-frame colors 

A variety of magnitude estimates are given for each galaxy in- 
cluded in the SDSS catalog. We chose the model magnitudes for 
deriving colors in our sample because they are calculated using 
best-fit parameters in r-band which are then applied to the other 
bands. Model magnitudes are therefore computed through the 
same aperture for all bands. In order to compute absolute mag- 
nitudes in each band the following corrections were applied: 

- A correction for Galactic dust extinction applying the 
reddening correctio ns computed by SDSS following 
ISchlegel et all (1 19981) . 

- The k-corr ection in each band was ca lculated using the code 
kcorrect dBlanton & Roweiil2007h . which determines the 
spectral energy distribution (SED) of the galaxy using SDSS 
photometry from a n onnegative linear combina tion of five 
templates based on the lBruzual & Charlotl d2003l) stellar evo- 
lution synthesis codes. 

- Absolute magnitudes were calculated using the updated dis- 
tances shown in Tab.Q] 

The internal extinction correctio n in B-band was calculated 
by IVerdes-Montenegro et al.l d2005 ) as a function of inc l inatio n 
and morphological type fo llowing de Vaucouleurs et all (1 199 lh . 
We used the Calzetti law (Calzet tiet al.ll2000h for deriving the 
SDSS-band internal extinctions based on the B-band correc- 
tions. We represented the rest-frame color (g-r) versus incli- 
nation before and after the extinction corrections and found that 
this correction was overestimated (the slope of the linear regres- 
sion changes from 0.17 to -0.21 for Sc gal axies and from 0.16 
to -0.13 for Sb galaxies, after corrections). Mast ers et al.l (2010) 
studied the effect of this correction on the colors through its 
dependency on inclination, spiral type and absolute magnitude. 
The best result for our (g-r) versus inclination plots is obtained 
using a linear combination of their relations as defined in equa- 
tion 3 in their paper (the slope of the (g-r) versus inclination re- 
lation after this correction is 0.025 for Sc and 0.023 for Sb galax- 
ies), while the expressions they give for each morphological type 
(bulge-disk ratio) produce an overcorrection in all our subsam- 
ples. We also note that the dependence of t he internal ext i nction 
correction on galaxy luminosity derived in Mas ters et alJ (2010) 
is not reproduced by our data. If AMIGA spirals are really more 
dynamically quiescent systems, it is perhaps not surprising that 
corrections derived from almost certainly less quiescent samples 
would yield different corrections. Moreover, this preliminary re- 
sult of lower extinction in AMIGA spirals can be understood as 
a sign of less dust due to our assumption that the star formation 
has b een lower in our isolated galaxies for all or most of their 
lives dLeon et alJl2008l) . 

Finally, we applied the correction of Mas ters et alJ d2010l) . 
based on their equation 3, to our Sa-Sc galaxies. We have not 
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Fig. 2. Distribution of Hubble types in the AMIGA-SDSS sam- 
ple. 
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Fig. 3. Distribution of the rest-frame (g-r) color as a function of 
Hubble type. The yellow points are the median values of each 
morphological type. The blue triangles are objects with asym- 
metry index IA=1 while the red open points represent the most 
asymmetric objects (IA=2) in our sample. The error bars repre- 
sent the median absolute deviation. 



corrected the morphological types earlier than (E/S0) and later 
than 6 (Scd/Sd/Irr) because their numbers are small and their 
(g-r) vs. inclination plots do not show any clear trend of color 
vs. i (the slope of the fits is ~ 0). 

5. (g-r) color as function of morphological type 

In Fig. [2] we presented the number of galaxies in the AMIGA- 
SDSS sample of each morphological subtype, and Fig.[3]shows 
the distribution of SDSS (g-r) colors for the AMIGA sample 
as a function of the morphological subtype. We calculated the 
median values rather than mean values, since they are less sen- 
sitive to outliers produced by an error in the color for one object 
or misclassification. We have used the median absolute devia- 
tion as a scatter measure since it represents a robust estimator of 
dispersion for median values. Not surprisingly the reddest me- 
dian values of (g-r) are found for the first four bins representing 



early-type galaxies although median (g-r) values remain essen- 
tially constant out to T=3 (Sb). Beginning with type Sb we see a 
decrease in median (g-r) as expected if this sequence reflects a 
uniformly decreasing contribution f rom an old stellar population 
(IStrateva et al.ll200UlBlanton et al.ll2003l) . 

Surprisingly the color distributions for our late-type spiral 
Sb-Sc population show some galaxies with colors even redder 
than the E/S0 subsample. Either AMIGA early-type galaxies 
are extraordinarily blue or these spirals have anomalous colors. 
Careful examination of these objects using the SDSS naviga- 
tion tool shows that colors of the reddest Sb-Sc galaxies are un- 
reliable and/or that the assigned galaxy types are incorrect. In 
the case of our Sc subsample (which is perhaps the easiest to 
classify) we find n=10 galaxies with (g— r)>0.7 (tipical colors 
of early-types). In 8 cases there is a problem with the color: 1) 
five are affected by a star projected on or very close to the galaxy 
(CIG 26 1 , 304, 649, 7 1 6, 1 1 0), 2) two (CIG 349, 946) belong to 
the 7 galaxies previously identified with a high degree of asym- 
metry (IA = 2), 3) one (CIG 988) is hig hly inclined and the ex- 
tinction correction (Mas ters et al.| [2010) could be maximally af- 
fected by uncertainty of the adopted inclination. Of the remain- 
ing 2 red Sc spirals, one (CIG 709) shows among the highest 
recession velocity in our sample (V>14000km/s). At this dis- 
tance the number of resolution elements in the galaxy images 
approaches that of POSS II. In the case of CIG 709, the image of 
SDSS shows a close object with unknown redshift which could 
be a minor companion causing some kind of asymmetry in the 
optical bands. Finally, CIG 807 may well be an Sc spiral but the 
blue disk is not resolved. The point of this exercise is to show 
that one can identify remarkably pure isolated (local) spiral sub- 
samples but only if the above pitfalls are taken into account. If 
these effects are not taken into account the average properties of 
the isolated subsamples can be blurred beyond statistical utility. 

Consideration of Sbc and Sb spirals reveal similar effects 
among the reddest galaxies that only detailed surface photometry 
could correct. Naturally similar problems can afflict the bluest 
galaxies in each spiral subtype. Colors for 2/4 of the bluest Sc 
spirals ((g-r)<0.4) are also likely affected by their nearby bright 
stars. These effects can be a source of serious scatter in color dis- 
tributions for any galaxy sample and illustrate the caution needed 
in using SDSS colors. On the other hand, while the objects with 
IA=2 are outside the normal trend of the median values, the color 
of the AMIGA galaxies with asymmetry index IA=1 agrees with 
the color of symmetric galaxies (IA=0). In the following analy- 
sis, objects with IA=2 have been removed but not objects with 
spurious colors, since they also exist in the samples that we have 
used as comparison. Nevertheless, we took account of their ef- 
fects through the error measures. 

We expect to find a smaller color dispersion for spiral sub- 
types in the AMIGA sample because these galaxies are mini- 
mally affected by environmental effects which apparently induce 
a larger color dispersion. This has been known since the first sta- 
tistical study of colors for galaxies in interact ing pairs and mul- 
tiplet s (i.e. they evolve more passively, e.g. Larso n & Tinslevl 
Il978l) . One of the goals of the AMIGA project is to provide a 
sample that could better characterize intrinsic galaxy properties 
and their dispersions. 

What is the main source of color dispersion for spiral sub- 
types if we have achieved our goal of minimizing effects of envi- 
ronmental nurture from our sample? Spirals contain a red bulge 
(redness depending on the nature of the bulge) and a blue disk 
making their global (g-r) color sensitive to the aperture used in 
estimating the g and r magnitudes. We checked the (g-r) color 
versus absolute magnitude relation for Sc galaxies in three sep- 
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arate redshift ranges. The results are presented in Fig. 3. We 
find a tendency for galaxies of fixed absolute magnitude to be 
(g-r)~0.08 redder at lower recession velocities implying that 
color measures in our sample could be increasingly affected by 
bulge light at lower recession velocities. This is a minor but sys- 
tematic effect. However, an aperture effect will also understimate 
the absolute magnitude of the galaxy, making it more difficult to 
estimate the amplitude of the color bias. 

Fig. |4] shows that the major source of color dispersion in 
Fig. [3] is connected with the color-luminosity trend. Sb-Sc 
color distributions show Gaussian distributions for each sub- 
type with FWHM=0. 15-0.2, spanning almost 2dex in luminos- 
ity (M r =-19 to -22.5). The effect is largest for Sc spirals in our 
sample where the color-luminosity trend appears to be steeper 
(Sc spirals with M r =-22.5 are ~0.25mag redder than Sc spirals 
with M, =-19). Fig. HI al so shows the (g- r ) vs. M , density dia- 
gram obtained from the iNair & Abraham! (|2()10) sample cover- 
ing the range 0.01<z<0.05 (SDSS DR8). This sample applied 
no morphology or environmental selection s so we see the bi- 
modality previously found for SDSS data ( Stra teva et alj feOOl ; 
Blan ton et al.l l2003). The distribution of our Sc subsample fol- 
lows the blue sequence with the recession velocity vs. color 
bias visible. The same bias in the col or as function of r edshif t 
has been found by using the sample of lNair &Abrahaml(l20Toh . 
Moreover, for a given luminosity and morphological type (Sc), 
the petrosian radius R9 given in the SDSS database for the 
INair & Abraham (2010) sample decreases as we go at lower red- 
shifts. The main trend of our data in the figure (and hence source 
of color dispersion) clearly involves the color-luminosity corre- 
lation. The overlap (i.e. green valley) rather than being a class 
of galaxies with intermediate colors may be entirely due to spi- 
rals with spurious red colors due to the effects discussed above. 
The (g-r) vs. M r diagram for a vetted sample of isolated galax- 
ies likely shows a much clearer dichotomy than other samples. 
There is no evidence for an astrophy sic ally significant green val- 
ley but we need more reliable colors to confirm this result. 

Fig-Hshows an equivalent diagram for the early-type (E/SO) 
part of the AMIGA sample. Despite their isolation they fall in 
the r ed sequence defined by early-types in richer environments 
(e.g. lNair & Abraham ( 2010)). They do not provide evidence for 
a green valley. Our early-type subsample shows no color trend 
with recession velocity. We notice only a slight trend between 
color and luminosity (galaxies with M r =-22.5 are ~0.1mag red- 
der than galaxies with M r =-19). In the case of early-types the 
FWHM of the Gaussian fit to the color distribution is ~0.14 
(similar to Sb galaxies color distribution). The bluest galaxies 
might be misc lassified as E/SO al though photometric study of a 
few of them (Marc um et al.1 12004) suggests that they may not be 
spirals. In addition, three of them (CIG 264, 981, and 1025) hav- 
ing HI emission were revised in order to check for the presence 
of peculiarities in the morphologies, such as dust lanes, optical 
shells, blue cores or st ar formation. A ll of them were confirmed 
as early-type galaxies (Espada 2006[). 



6. Color dependence on environment 

Some studies have shown that galaxies of a fixed morphol- 
ogy in higher dens ity environments are redder (iBall et alJl2008l 
Skibba et al. 2009) and the reason might be that dense environ- 
ments suppress star formation (Balogh et al. 1998). On the other 
hand isolated galaxies are likely to show rather passive star for- 
mation leading one to expect spirals in richer environments (e.g. 
pairs, even the isolated pairs usually included in field galaxies 
samples) to show more active star formation and hence bluer 
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Fig. 4. (g-r) color-magnitude diagram for the Sc galaxies in the 
AMIGA sample. The blue diamonds are objects at z<0.02, the 
yellow triangles are galaxies at redshift 0.02<z<0.04, and the 
pink squares are those objects at z>0.04. Th e grey-scale rep- 
resent s the density diagram obtained from the Nair & Abraham 
(I2010I) sample at 0.01<z<0.05, using the data of DR8. 
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Fig. 5. Same as Fig.|4]but for the E/SO galaxies in the AMIGA 
sample. 



colors dPatton et al.l201 lh . lLarson & Tinslevl(ll978l) found larger 
color dispersion in an extreme nurture sample, whose evolution 
is completely dominated by external effects. The key word here 
is dispersion: colors of galaxies in pairs and groups are not sys- 
tematically bluer but they show a larger color dispersion. The 
AMIGA sample of galaxies involves some of the most isolated 
objects in the local Universe and are minimally affected by ex- 
ternal processes, at least during the past few Gyr. We therefore 
expect that AMIGA galaxies of a given morphological type will 
show minimal color dispersion and the above results support that 
expectation. 
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Table 2. Median (g-r) colors as function of morphological type. The errors represent the median absolute deviation. 



Type 


T 


AMIGA 


NAIR 


EFIGI 


CPG 


CPG 












(WID) 


(CLO) 


E 


-5 


0.72+0.06 


0.78±0.03 


0.78+0.03 


0.79+0.03 


0.76+0.03 


E 


-4 


0.77±0.02 


- 


0.78±0.02 


0.80+0.08 


0.79±0.04 


E/SO 


-3 


0.75±0.04 


0.76±0.05 


0.77±0.04 


0.79±0.06 


0.77+0.07 


SO 


-2 


0.75+0.04 


0.76±0.04 


0.76±0.04 


0.78±0.06 


0.77+0.06 


SO 


-1 


- 


- 


0.78±0.06 


0.72±0.05 


0.73±0.09 


SO/a 





0.74±0.07 


- 


0.76±0.07 


0.77+0.05 


0.78±0.04 


Sa 


1 


0.77±0.05 


0.71+0.06 


0.73±0.05 


0.72±0.11 


0.71 ±0.09 


Sab 


2 


0.74±0.05 


0.69±0.07 


0.72±0.07 


0.71+0.10 


0.67±0.15 


Sb 


3 


0.71+0.06 


0.67±0.08 


0.71+0.08 


0.71±0.13 


0.69±0.12 


Sbc 


4 


0.65±0.09 


0.61±0.08 


0.66±0.07 


0.63±0.12 


0.59±0.14 


Sc 


5 


0.57±0.08 


0.56±0.08 


0.62±0.09 


0.69±0.12 


0.51±0.15 


Scd 


6 


0.49±0.06 


0.46±0.07 


0.58±0.09 


0.55±0.11 


0.51±0.17 


Sd 


7 


0.42+0.06 


0.42±0.06 


0.47±0.08 


0.34±0.17 


0.43±0.10 


Sdm 


8 


0.35±0.05 


0.41±0.07 


0.44±0.12 


0.48±0.07 


0.30±0.12 


Sm 


9 




0.36±0.09 


0.40±0.16 




0.56±0.13 


Im 


10 


0.24±0.03 


0.33±0.10 


0.29±0.09 




0.29±0.12 



Table [2] provides a quantitative comparison of (g-r) me- 
dian colors in our sample and in samples involving denser en- 
vironments. We used t hree catalogs for this comparison: 1) 
iNair & Abraham (2010) which includes a detailed visual classi- 
fications for 14,034 galaxies in the SDSS DR4. We selected only 
objects with morphological classification available and redshift 
0.01<z<0.05 (8976) to better matc h our AMIGA-SPS S sam- 
ple ( 98%); 2) the EFIGI catalog dBaillard et al.ll201ll) which 
provides detailed morphological information for a sample of 
4458 PGC galaxies also using SDSS DR4. Both samples include 
galaxies in a wide range of environments; 3) We also compare 
with t he catalog of isolated pairs of galaxies (CPG. iKarachentsevI 
|1972|) . composed by 1206 objects. It was visually compiled, like 
CIG, using POSS II and applying an isolation criterion. 

We used SDSS DR8 photometric data for all three sam- 
ples and derived rest-frame colors in the same way as for 
our sample. For the CPG sample, there are 916 galaxies with 
SDSS photometry. We calculated the projected distance between 
galaxies i n each pair using red shifts taken from HyperLeda. 
Following IXu & Sulentid (Il99ll) we separated this sample into 
close-interacting (CLO - galaxies with separation SEP< 2Mpc) 
and wide (WID) pairs (SEP>2Mpc). We checked the distribution 
of redshift and absolute magnitudes in the r-band for each sam- 
ple and found good agreement with the whole AMIGA-SDSS 
sample. In the cases of EFIGI and CPG catalogs, we removed 
the objects with recession velocities lower than 1500 km/s as we 
did for our AMIGA sample, due to the determination of the iso- 
lation parameters dVerlev et al.ll200 7b). We have also removed 
objects fainter than M r >-17 since we do not have counterparts 
in the AMIGA-SDSS sample. The final samples of compari- 
son have 8879 (Nair), 3294 (EFIGI), and 839 (CPG) galaxies. 
In Fig. [6] we present the absolute magnitude distribution in the 
r-band for each sample, as well as the probability given by the 
Kolmogorov-Smirnov two sample test that these distributions 
are indistinguishable to the AMIGA-SDSS one. This probabil- 
ity is larger than 90% for the CPG and Nair samples, while the 
lowest value is found in the case of the EFIGI sample (63%). 

We found good agreement for the colors of early-type galax- 
ies in all samples. However, in all cases (from T=-5 to T=0), 
the median values for the AMIGA sample are a little bluer 
(but within the errors) than the other samples. Because the 
color-luminosity relation has a trend (brighter galaxies are red- 
der), this bluer color could be presumably due to the fact that 
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Kolmogorov-Smirnov test: 

Prob AMIGA-CPG = 0.999 
Prob AMIGA— Noir= 0.91 2 
Prob AM IGA— Ef ig I = 0.633 
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Fig. 6. Distribution of absolute magnitudes i n the r-band for 
the wh ole AMIGA-SDSS sample, and for fhe lNair & Abraham 
(2010), Efigi, and CPG samples described in the text. 



the E/S0 population of A MIGA is fainter than any other sam- 
pleJjSulentic et al. 2006). The median values of colors for the 
INair & Abrahaml ( 1201 Oft and close pairs samples, for the (Sb-Sc) 
spirals, are consistent but sli ghtly bluer than our sam ple. In the 
case of the comparison wifh lNair & Abraham! d2010t) . the differ- 
ences in color may be produced by the different morphological 
classifications used in both samples. We find their morphologies 
to be earlier than ours, with a mean deviation of ~ 1 .5 for each 
Hubble type. In order to test this possibility, we selected a sub- 
sample composed by the 1 42 common objects between AMIGA 
and the lNair & Abraham] d2010h sample, and calculated the me- 
dian values of (g-r) obtained using each morphological classifi- 
cation. We found that the difference between both median colors 
in the common sample was consistent with that found using the 
full samples. 

Nevertheless, the differences in color of close pairs seem to 
be more robust because the colors of wide pairs are as red as 
the colors of isolated galaxie s. The scatter is also larg er for the 
close pairs, in agreement wifh lLarson &Tinslevldl978h .InFig.l7l 
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Fig. 7. Distribution of absolute magnitudes in the g-band (left), r-band (centre) and (g-r) color (right) for all Sb (top), Sbc (middle) 
and Sc (bottom) galaxies. The solid lines represent the AMIGA-SDSS sample and the dashed lines are the distribution of the CPG 
close pairs. 



we present the distribution of absolute magnitudes in r and g- 
bands as well as (g-r) color for three morphological types (Sb, 
Sbc and Sc) of the AMIGA-SDSS sample. We have also rep- 
resented the same distributions for the close pairs of the CPG 
sample. While the color distribution of each morphological type 
has gaussian distribution for our sample, the color histogram for 
the CPG sample follows a non gaussian distribution. The loss 
of Gaussianity in the distribution of color can be interpreted as 
an effect of the environmental nurture that occurs in a sample of 
pairs. 



We have also investigated the change in the median colors 
when we rejected objects with asymmetries (IA=1) and the dif- 
ferences are negligible. 

7. Conclusions 

We have taken a first look at the colors of isolated galaxies using 
SDSS g and r magnitudes. A first look is also inevitably an intro- 
duction into the pitfalls of using the SDSS automated measures 
of magnitudes. We find that the color distributions of morpho- 
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logical subtypes can be well described as Gaussian distributions 
with FWHM (g-r)=0. 1-0.2. This is expected for samples where 
effects of environmental nurture have been minimized, and is 
supported by the fact that this Gaussianity was not observed in 
the sample of galaxy pairs. The majority of the color dispersion 
for spiral subtypes is caused by a color-luminosity correlation, 
with more massive spirals showing redder colors. In fact most 
of the Sb spirals in our sample concentrate at the upper region 
of the blue cloud, below the red sequence. We see little evidence 
for a green valley in our sample with most spirals redder than 
(g-r)=0.7 having spurious colors (~80%). We find a preliminary 
difference in the color of the isolated and paired spiral galax- 
ies. The median value of (g— r) seems to be bluer when inter- 
actions come into play. Nevertheless, no difference is found in 
the (g-r) color of early-type galaxies, which suggests that the 
bluer color of spirals in pairs is presumably due to interaction 
enhanced star formation. Our sample of isolated galaxies gives 
a median absolute deviation in color which is lower than that in 
pairs of galaxies, where a more active star formation and, per- 
haps, a higher dust diffusion caused by the interaction, are also 
sources of color dispersion. The mean colors and dispersions for 
isolated galaxy subtypes (especially E/SO and Sb-Sc) are likely 
the best nurture-free measures so far obtained. 
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